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Introduction 

An important element in the National Aero-Space 
Plane (NASP) research program at the Langley 
Research Center is the development of the analytical 
tools or models that are needed for analyzing and 
gaining insight into the combustion problems associ- 
ated with propulsion systems required to operate at 
flight speeds up to Mach 25. Recently, special at- 
tention has focused on the role that chemical kinetic 
effects can have on the performance of high Mach 
number combustion systems. Mechanisms describing 
the detailed chemical kinetics of hydrogen oxidation 
have been developed by a number of authors, and 
many of these mechanisms have been reviewed by 
Dougherty and Rabitz (1980). In the present paper a 
detailed chemical kinetic mechanism for the combus- 
tion of hydrogen is presented and discussed. The pur- 
pose of this study was to assemble a chemical mecha- 
nism that reproduced, with reasonable accuracy, the 
observed behavior of the hydrogen-air system as de- 
termined in shock-tube studies and in laminar flame 
studies. A second objective was to examine the effect 
of chemical kinetics on combustion at representative 
scramjet combustor conditions. Particular attention 
was given to the effect of uncertainties in the rate co- 
efficients on the calculated ignition delay times, burn- 
ing velocities, and combustor performance. 

Symbols 

A/A 0 area ratio along combustor 

a constant in mixing schedule equation 

D self-diffusion coefficient for molecular 

oxygen 

Dij binary diffusion coefficient for species i 
and j 

Fi,Fj diffusion factors for species i and j 

kj rate coefficient for reaction j 

M flight Mach number 

MWj molecular weight of species i 

p pressure, atm 

R universal gas constant 

T absolute temperature, K 

V initial velocity at combustor entrance 

r ) m ix mixing efficiency 

r/ c *. chemical kinetic efficiency 

<f> fuel-air equivalence ratio 


Hydrogen-Air Reaction Mechanism 

In the assembly of the reaction mechanism, re- 
actions were included that involve all the impor- 
tant species in the hydrogen-oxygen system — H 2 , 
0 2 , H, O, OH, H 2 0, H0 2 , and H 2 0 2 . In ad- 
dition, reactions were included for the species N, 
NO, and HNO, which may become important in 
the hydrogen-oxygen-nitrogen system at the high 
Mach number conditions ( M > 12). The proposed 
hydrogen-air combustion mechanism and rate coeffi- 
cients are given in table I. The reactions in the mech- 
anism are based for the most part on the reactions 
and rate coefficients recommended by Baulch et al. 
(-1972), Baulch, Drysdale, and Horne (1973), Dixon- 
Lewis and Williams (1977), Dixon-Lewis (1979), and 
Dougherty and Rabitz (1980). The third-body effi- 
ciencies were taken from Gay and Pratt (1971) and 
Dixon-Lewis (1979). The rate coefficient expressions 
assigned to reactions (2) and (9) represent values, 
optimized from this study, that gave the best overall 
agreement between calculated kinetic results and ex- 
perimental data. The refinement of the hydrogen-air 
reaction mechanism is discussed in the next section. 

Mechanism Evaluation and Refinement 

The hydrogen-air reaction mechanism was refined 
and evaluated through comparison of calculated ki- 
netic results with experimental data. The refinement 
involved the adjustment of the rate coefficients for 
certain reactions to obtain the best agreement with 
experimental data. The experimental data selected 
for the comparison were the shock-tube ignition de- 
lay data reported by Slack (1977), the laminar burn- 
ing velocity data assembled by Warnatz (1981), and 
the laminar burning velocity data reported by Milton 
and Keck (1984). These data were selected for the 
comparison because they were obtained from experi- 
ments with “real” hydrogen-air mixtures as opposed 
to diluted or simulated hydrogen-air (e.g., H 2 - 02 -Ar) 
mixtures. 

Comparison With Shock-Tube Data 

Slack (1977) measured ignition delay times be- 
hind the reflected shock for stoichiometric hydrogen- 
air mixtures near the second explosion limit and 
for stoichiometric mixtures containing small amounts 
of nitric oxide. The data selected for the com- 
parison included the ignition delay data at pres- 
sures of 0.5, 1.0, and 2.0 atm for stoichiometric 
hydrogen-air mixtures and the ignition delay data at 
2.0 atm for stoichiometric hydrogen-air mixtures con- 
taining 0.50 and 2.25 percent nitric oxide. To com- 
pare the kinetic behavior predicted by the reaction 
mechanism with the experimental results, numerical 



simulation of the shock-tube experiments was car- 
ried out with the computer code described by McLain 
and Rao (1976). To simulate the conditions behind 
the reflected shock wave, the code was operated in 
a constant-volume mode. The ignition delay time 
was defined as it is in the shock-tube experiments; 
that is, it was defined to be the elapsed time be- 
tween the heating of the gas mixture by the reflected 
shock wave and the sudden pressure increase due to 
combustion. 

The rate coefficients for the reverse reactions 
in the mechanism were calculated within the com- 
puter code with the forward rate coefficients given in 
table I and the appropriate thermochemical data. 
The thermochemical data for the hydrogen, oxy- 
gen, and nitrogen species were taken from the 
JANAF (1971) tables. 

The calculated ignition delay times were very 
sensitive to the rate coefficients assigned to reactions 
(2) and (9), 

(2) H + O 2 — ► OH + O 
(9) H + 0 2 + M -> H0 2 + M 

particularly for the 2.0-atm condition. The sensi- 
tivity of the calculated ignition delay times to reac- 
tions (2) and (9) is illustrated by the results shown 
in figure 1. The solid curve in this figure was com- 
puted with the rate coefficients listed in table I, 
and the dashed curves were computed with the al- 
tered rate coefficients noted on the figure. Cal- 
culated burning velocities (discussed subsequently) 
were also very sensitive to these rate coefficients, 
particularly for some rich mixtures. The rate co- 
efficients initially assigned to reactions (2) and (9) 
were taken from Baulch et al. (1972). Baulch rec- 
ommended == 2.2 x 10 14 exp(— 16800/RT) with 
an uncertainty of ±50 percent for the temperature 
range 300 to 2000 K. The best overall agreement be- 
tween the calculated and experimental results was 
obtained with the rate coefficients for reactions (2) 
and (9) given in table I. The adjusted rate coeffi- 
cient for reaction (2) was 18 percent larger than the 
recommended rate coefficient, and the adjusted rate 
coefficient for reaction (9) was 5 percent larger than 
the recommended rate coefficient. The adjusted ex- 
pressions were within the accepted uncertainty for 
the reactions. 

The results of the computer simulations of the 
shock-tube experiments for stoichiometric hydrogen- 
air mixtures at 0.5, 1.0, and 2.0 atm are given in 
figures 2(a) to 2(c) together with the corresponding 
experimental results. The calculated ignition delay 
times were within 20 percent of the experimental 
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Figure 1. Sensitivity of calculated ignition delay 

times to rate coefficients for reactions (2) and (9). 

results over the entire range of temperatures and 
pressures examined. 

At high flight Mach numbers (M > 12), reac- 
tions involving nitric oxide are important. Slack and 
Grillo (1977) investigated the effect of nitric oxide ad- 
dition on the ignition of stoichiometric hydrogen-air 
mixtures and determined that the presence of small 
amounts of nitric oxide (less than 5 percent) reduced 
the ignition delay times. The high-temperature 
chemistry of the nitrogen-oxygen system is repre- 
sented in the mechanism by reactions (20) to (32). 
The experimental results of the study at 2.0 atm for 
the addition of 0.5- and 2.25-percent nitric oxide are 
given in figure 3 together with the results calculated 
with the proposed reaction mechanism. The agree- 
ment between the calculated and experimental re- 
sults is good. The sensitizing effect of the nitric oxide 
is primarily due to 

(30) H0 2 + NO -► N0 2 + OH 

(31) H + NO 2 — NO + OH 

which convert the HO 2 , a chain- terminating species 
during the ignition delay period, to the very reac- 
tive OH radical. The rate coefficients assigned to 
reactions (22) to (33) were based primarily on the 
expressions recommended by Slack (1977) and ad- 
justments were not required. 

Comparison With Burning Velocity Data 

The ability of a reaction mechanism to adequately 
predict burning velocities is an important, and gen- 
erally more comprehensive, test of the mechanism at 
conditions more representative of a combustion en- 
vironment. The data selected for the comparison 
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(b) p = 1.0 atm. 



Figure 2. Calculated ignition delay times compared with experimental results for stoichiometric H 2 
mixtures. 
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Figure 3. Calculated ignition delay times compared with experimental results for stoichiometric H 2 -air 
mixtures containing nitric oxide, p = 2.0 atm. 


were the burning velocities assembled by Warnatz 
(1981) for hydrogen-air flames at pressures of 1 atm 
and an initial temperature of 298 K for equivalence 
ratios between 0.5 and 6.0. The burning velocity 
data reported by Milton and Keck (1984) for stoi- 
chiometric hydrogen-air mixtures at temperatures up 
to 515 K and pressures up to 7 atm were also used. 
Burning velocities were calculated with the premixed 
one-dimensional flame (PROF) code developed origi- 
nally by Kendall and Kelly (1978) for premixed, lam- 
inar, laterally unconfined flames. The PROF code 
models the axial diffusion of heat and species, as well 
as the chemical kinetic processes, that occur within 
a one-dimensional flame. The code treats each reac- 
tion as reversible and calculates the rate coefficient 
for the reverse reaction using the forward rate coeffi- 
cient and the appropriate equilibrium constants. The 
thermochemical data were obtained from the JANAF 
(1971) tables. 

The governing equations in the flame model are 
developed through integration of the steady-state 
two-dimensional species, momentum, and energy 
equations across a plane perpendicular to the axis. 
This results in a set of one-dimensional flame equa- 
tions in terms of bulk gas properties. Species trans- 
port equations in the code use binary diffusion coeffi- 
cients D{j, which are approximated by the empirical 
relation 


D 

Di i ~ p-F- 

r t r 3 

where D is a reference self-diffusion coefficient and 
Fi and Fj are diffusion factors for species i and 
j. Pressure and temperature dependence of D{j are 
incorported into D so that F t and Fj are independent 
of temperature and pressure. The reference species 
is molecular oxygen, for which 

D = 1.7 x 10- 5 r 1 - 65 V 1 cm 2 - sec -1 

where the pressure and temperature are in atmo- 
spheres and kelvins. The diffusion factor for species 
i is determined from the empirically derived relation 

Fi = (MWj/26) 0461 

where MWj is the molecular weight of species i. 
Additional details on the development of the binary 
diffusion coefficient relation are given in Kendall 
and Kelly (1978) and Bartlett, Kendall, and Rindal 
(1968). 

As noted previously, burning velocity data were 
also used to evaluate and optimize the proposed 
hydrogen-air reaction mechanism. The data used for 
this purpose were those compiled by Warnatz (1981) 
for various hydrogen-air mixtures for a pressure of 
1.0 atm and an unburned gas mixture temperature 
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of 298 K. These data are shown in figure 4 together 
with curves that illustrate the sensitivity of the cal- 
culated burning velocities to various reactions. The 
initial mechanism (without any adjustments to the 
rate coefficients for reactions (2) and (9)) correctly 
predicted the position of the maximum burning ve- 
locity (H 2 mole fraction ss 0.4) but underpredicted 
the actual value by almost 20 percent. Sensitivity 
studies with the mechanism showed that the value of 
the maximum burning velocity was most sensitive to 
the rates of reactions (2) and (9) and, to a lesser ex- 
tent, to the net rate of the H02-consuming reactions 
(10) to (15): 

(2) H + 0 2 -> OH + O 

(9) H + 0 2 + M — *• H0 2 + M 

(10) HO 2 + H-+H 2 + O 2 

(11) H0 2 +H^0H + 0H 

(12) H0 2 + H^H 2 0 + 0 

(13) H0 2 + O -f 0 2 + OH 

(14) H0 2 + OH — ► H 2 0 + 0 2 

(15) H0 2 + H0 2 -+ H 2 0 2 + 0 2 

The importance of the chemistry of the H0 2 radi- 
cal in the reaction mechanism is illustrated in fig- 
ure 4. When the rate coefficient for reaction (9) was 
decreased by a factor of 10, effectively cutting off 
the production of the H0 2 radical, the calculated 
burning velocities were significantly smaller than the 
experimental data. In the hydrogen-air flame the 
HO 2 radical is a very important species in promoting 
flame propagation through reaction (11), which gen- 
erates the hydroxyl radical OH. The hydroxyl radi- 
cal then reacts with the molecular hydrogen through 
reaction (4): 

(11) HO 2 + H-OH + OH 
(4) OH + H 2 ^H 2 0 + H 

The results of the sensitivity studies with the sim- 
ulated shock-tube and flame experiments indicated 
that the best agreement between the calculated re- 
sults and experimental data could be achieved by a 
simultaneous increase in the rate coefficients for re- 
actions (2) and (9) while the rate coefficients for re- 
actions (10) to (15) were kept fixed. Adjustments 
to the rate coefficients were kept within the reported 
uncertainty limits. The solid curve shown in figure 4 
represents the burning velocities calculated with the 
mechanism listed in table I. The mechanism correctly 
predicted the position of the maximum burning ve- 
locity and the calculated values were within 5 percent 
of the experimental data. 
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Figure 4. Calculated burning velocities compared 

with experimental data for H 2 -air mixtures at 

p = 1.0 atm and unburned gas temperature of 
298 K. 

The results shown in figures 5 to 7 are compar- 
isons of calculated burning velocities with the data 
reported by Milton and Keck (1984). Milton and 
Keck measured laminar burning velocities for stoi- 
chiometric hydrogen-air mixtures using a spherical 
constant-volume combustion bomb which gave re- 
sults over a range of unburned gas temperatures and 
pressures. The results shown in figure 5 illustrate 
the effect of the initial unburned gas temperature on 
the burning velocity for a stoichiometric mixture at 
p = 1.0 atm. The calculated results were in good 
agreement at 300 K and about 20 percent higher at 
400 K. The results shown in figure 6 illustrate the 
effect of pressure on the burning velocity at an un- 
burned gas temperature of 298 K. The calculated re- 
sults were within 8 percent of the experimental data. 

Most of the hydrogen-air burning velocity data 
generated by Milton and Keck (1984) were for a range 
of unburned gas temperatures and pressures. These 
results fall within the data band shown in figure 7. 
The calculated burning velocities were in good agree- 
ment with the experimental results for pressures up 
to about 5.0 atm. At higher pressures the experimen- 
tal burning velocities were up to 15 percent larger 
than the calculated results. Attempts to provide 
better agreement at the higher pressures (by vary- 
ing rate coefficients) without affecting the results at 
other conditions were not successful. No explanation 
can be offered for the difference between the calcu- 
lated and experimental results at the high pressures 
and temperatures. 

Application to Scramjet Combustion 

Calculations have been carried out to examine the 
effect of chemical kinetics on combustion at repre- 
sentative scramjet combustor conditions. The com- 
bustor conditions that were considered correspond 
to flight Mach numbers of 8, 16, and 25. The math- 
ematical model used in this study was based on a 
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Figure 5. Comparison of calculated burning veloci- 
ties with experimental results for p = 1.0 atm. 
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Figure 6. Comparison of calculated burning veloci- 
ties with experimental results for T = 298 K. 
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Figure 7. Comparison of calculated burning veloci- 
ties with experimental results for range of tem- 
peratures and pressures. 


one-dimensional flow reactor which solved the mass, 
momentum, and energy conservation equations for 
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a steady flow system. The model was based on 
the computer codes described by Bittker and Scullin 
(1972) and McLain and Rao (1976) and was capable 
of treating reactions in a flowing system with variable 
geometry. To simulate the combustion process it was 
assumed that at the entrance to the combustor the 
hydrogen, air, and any ignition source constituents 
are in thermal equilibrium but not fully mixed on 
a microscale. A mixing routine was included in the 
model, which allowed the unburned fuel and air to be 
mixed at a prescribed schedule along the combustor. 
The fuel and air were mixed; they were then per- 
mitted to react. The mixing schedule was fixed to 
maintain a stoichiometric reaction zone until either 
the fuel (when 0 < 1) or the air (when 0 > 1) was 
fully mixed. 

The conditions at the entrance to the combustor 
are given in table II for the three flight Mach num- 
bers examined. Column T is the temperature of the 
hydrogen-air mixture and V is the initial velocity of 
the mixture. These initial conditions are for a com- 
bustor in which the fuel is injected perpendicular to 
the direction of the air flow. The combustor geome- 
try associated with each flight Mach number is shown 
in figure 8 in terms of the combustor area distribu- 
tion. Each combustor configuration had two sections 
with different area distributions. The first section 
was 50 cm long and the second section was 100 cm 
long. Even though the configurations in figure 8 have 
been referred to as a combustor, subsequent discus- 
sions in this report refer to the first section as the 
combustor section and the second section as the noz- 
zle section. 

The assumed mixing schedules for each condition 
are shown in figure 9. The mixing schedule was 
input into the computer code in the form rj m \x = 
1 — exp(— ax), where a is a constant, x is the distance 
in centimeters, and rj m ^ is the mixing efficiency. 
The mixing schedule was obtained from a correlation 
developed by Anderson and Rogers (1971) and by 
Anderson (1974). The correlation, which is a func- 
tion of fuel-air equivalence ratio and fuel injector 
spacing, was based on experimental results from 
studies of cold flow mixing. When the overall equiv- 
alence ratio is less than or equal to one, the mixing 
efficiency ? 7 m i x is defined as the fraction of total hy- 
drogen that has been mixed with a stoichiometric 
amount of air and allowed to react. When the equiv- 
alence ratio is greater than one, the mixing efficiency 
is defined as the fraction of total air that has been 
mixed with a stoichiometric amount of hydrogen. 

The effect of chemical kinetic factors on the com- 
bustion process was evaluated by comparison of the 
calculated specific internal thrust for the case in 
which the combustion chemistry was described by 
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Figure 8. Combustor area distributions. 

nonequilibrium chemistry with that for the case in 
which the combustion chemistry was described by 
equilibrium chemistry. The specific internal thrust 
is defined as the internal thrust divided by the total 
mass flow of fuel and air. The results of the study are 
expressed in terms of the ratio of the internal thrust 
for the nonequilibrium case to the internal thrust for 
the equilibrium case. In this report this ratio is called 
the chemical kinetic efficiency 

The results of the initial studies for the Mach 8 
condition revealed that self-ignition of the hydrogen- 
air mixture did not occur and that an ignition source 
was needed. Consequently a study was carried out 
to determine what effect the quantity of the ignition 
source had on the combustion process. The results of 
this study are shown in figure 10, in which the chem- 
ical kinetic efficiency is plotted against the ignition 
source fraction. When the equivalence ratio is less 
than or equal to one, the ignition source fraction is 
defined as the fraction of the total moles of hydrogen 
that is allowed to react to equilibrium with a stoi- 
chiometric amount of air and to act as an ignition 
source for the remaining hydrogen and air. When 
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Figure 9. Mixing schedules for each Mach number. 

the equivalence ratio is greater than one, the igni- 
tion source fraction is defined as the fraction of the 
total amount of air that is allowed to react to equilib- 
rium with a stoichiometric amount of hydrogen. The 
purpose of the ignition source is to simulate a flame- 
holding region that develops around the fuel injector 
or behind a step in the combustor. From a modelling 
perspective the ignition source material was assumed 
to adiabatically mix with remaining fuel and air, a 
condition which results in an initial mixture with a 
higher temperature and the presence of free radicals. 
The mixing schedules given in figure 9 are used for 
the remaining fuel and air. 

The results in figure 10 show that for the Mach 
8 condition an ignition source fraction of about 0.1 
was required to achieve a reasonable chemical kinetic 
efficiency near 0.9. For the Mach 16 and Mach 25 
conditions the presence of an ignition source did 
not result in a significant increase in the chemical 
kinetic efficiency. For these conditions the initial 
temperature and pressure were high enough to allow 
self-ignition to occur. For the Mach 25 condition 
sufficient oxygen atoms were present in the air and 
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Figure 10. Effect of ignition source fraction on chemical kinetic efficiency. 


the additional free radicals provided by the ignition 
source had a very small effect on the chemistry. 

It is evident from the results shown in figure 10 
that chemical kinetic effects were present because the 
chemical kinetic efficiency was less than 1.0 at the 
conditions examined. To determine which reaction 
processes were responsible for the kinetic limitations, 
a parametric study was performed in which the rate 
coefficients for certain reactions were varied within 
specified limits. Initially, we varied the rates of var- 
ious groups of reactions rather than the rates of in- 
dividual reactions. The groups that were considered 
included reactions (2) to (5), 

(2) H + O 2 — ► OH + O 

(3) O + H 2 -* OH + H 

(4) OH + H 2 H 2 0 + H 

(5) OH + OH-+H 2 O + O 

which are the primary chain branching and propa- 
gating reactions; reactions (6) to (8) and (20), 


(6) H + OH + M — ► H 2 0 4- M 

(7) H + H + M-+H 2 + M 

(8) H + O + M — ► OH + M 
(20) O + O + M-+O 2 + M 

which are the primary recombination reactions; and 
reaction (9) 

(9) H + 0 2 + M -> H0 2 + M 

which is the primary source of HO 2 in the hydrogen- 
air reaction mechanism. Reactions involving HO 2 
were important in both the ignition delay calcula- 
tions and the burning velocity calculations. The 
parametric studies were carried out for hydrogen-air 
mixtures with an ignition source fraction of 0.10. The 
results of the studies are shown in figures 11 to 13. 
In the figures the temperature distributions along the 
combustor are plotted for the various cases. Results 
are given for cases in which the rate coefficients in 
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Figure 11. Chemical kinetic effects at Mach 8 condition for <j> = 1 and ignition source fraction of 0.10. 



Figure 12. Chemical kinetic effects at Mach 16 condition for (j) = 1 and ignition source of 0.10. 
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Figure 13. Chemical kinetic effects at Mach 25 condition for 4 ) = 2 and ignition source of 0.10. 


table I for reactions (2) to (5) were varied by a fac- 
tor of 2, those for the recombination reactions (6) to 
(8) and (20) were varied by a factor of 5, and that 
for reaction (9) was varied by a factor of 10. Also 
shown in each figure is the chemical kinetic efficiency 
for each case. The reference case in each figure refers 
to the result obtained with the rate coefficients given 
in table I. 

The results for the Mach 8 condition are shown 
in figure 11. The results show that all the reactions 
examined were important and contributed to the ki- 
netic effect. It is evident that at the Mach 8 condition 
there was an ignition delay which was controlled pri- 
marily by the rates of reactions (2) to (5). A decrease 
in the rate of these reactions produced a significant 
decrease in the chemical kinetic efficiency, whereas 
an increase in the reaction rates produced only a 
slight increase in the chemical kinetic efficiency. Vari- 
ation of the rates of the recombination reactions (6) 
to (8) and (20) also influenced the chemical kinetic 
efficiency. The results of additional parametric stud- 
ies revealed that the recombination reactions had 
the largest effect in the combustor section, with 
reaction (6) 

(6) H + OH + M H 2 0 + M 

having the largest influence. Variation of the recom- 
bination rates in the nozzle section did not produce 
significant changes in the chemical kinetic efficiency 
even though some reaction continued in this section. 


A significant result of the Mach 8 study was 
the demonstration that the H0 2 chemistry was an 
important element in the reaction mechanism. A 
reduction in the rate of reaction (9) by a factor of 10 
(see curve 7 of fig. 11) produced a 15-percent decrease 
in the chemical kinetic efficiency, while an increase 
in the rate by a factor of 10 produced a 5-percent 
increase in the chemical kinetic efficiency. This effect 
can be explained in terms of reactions (9) and (11): 

(9) H + 0 2 + M — ► H0 2 + M 
(11) H + H0 2 -+OH + OH 

Reaction (9) is the primary source of the H0 2 radi- 
cal and reaction (11) is the primary chain branching 
reaction involving the H0 2 radical. Reducing the 
rate of reaction (9) results in less H0 2 production 
and, consequently, less chain branching. In addition, 
reactions (9) and (11) are exothermic; therefore, a re- 
duction in reaction (9) results in less heat release and 
consequently a lower temperature, which affects the 
rate of the important endothermic chain branching 
reaction 

(2) H + 0 2 OH + O 

When the rate of reaction (9) is increased, the effect 
on the chemical kinetic efficiency is not as large as 
the effect when the rate of reaction (9) is decreased 
because even though more H0 2 radical is produced, 
the conversion of H0 2 to OH is limited by the rate 
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of reaction (11) and competition from reactions (10) 
and (14) 

(10) H0 2 + H H 2 + 0 2 

(14) H0 2 + OH -► H 2 0 + O 1 2 

For the Mach 8 condition the H0 2 chemistry is 
important primarily in the combustor section and is 
less important in the nozzle section, even though the 
reactions continue. 

The results for the Mach 16 condition are shown 
in figure 12. For this condition, varying the rates of 
reactions (2) to (5) did not affect the chemical kinetic 
efficiency, and this lack of effect indicates the over- 
all combustion process was controlled by other reac- 
tions. The recombination reactions (6) to (8) and 
reaction (9) had the largest influence on the combus- 
tion process and this effect primarily occurred in the 
combustor section. In the nozzle section the kinet- 
ics of the three-body reactions were essentially frozen 
while the two-body reactions continued. The results 
given in figure 12 again illustrate the importance of 
the H0 2 chemistry and the need to include the H0 2 
chemistry in the mechanism. 

The results for the Mach 25 condition are shown 
in figure 13. Varying the rates of reactions (2) to 
(5) did not affect the chemical kinetic efficiency. The 
results of additional studies revealed that the con- 
trolling reactions in the combustor were the recom- 
bination reactions (6) to (8), with reaction (6) being 
the most important in the combustor section. In the 
nozzle section the recombination reactions were vir- 
tually frozen. 

Discussion of Results 

The results of the study suggest that chemi- 
cal kinetic effects can be important at representa- 
tive scramjet combustor conditions and that combus- 
tor models which use nonequilibrium chemistry are 
preferable to models that assume equilibrium chem- 
istry. For the conditions examined the; results also 
show the importance of including the H0 2 chemistry 
in the mechanism. Results obtained with a simplified 
model that does not include the H0 2 chemistry can 
be significantly different from the results obtained 
with a model that includes H0 2 chemistry. This 
difference is illustrated in figures 14 and 15, which 
show the temperature profile and the chemical ki- 
netic efficiency for the Mach 8 and Mach 16 condi- 
tions. The solid curves are the results based on the 
reaction mechanism and the rate coefficients listed 
in table I. For the Mach 8 condition omitting the 
H0 2 and H 2 0 2 reactions from the model resulted in 
a 29-percent reduction in the chemical kinetic effi- 
ciency, while for the Mach 16 condition there was a 


6-percent reduction. The reactions primarily respon- 
sible for this effect were reactions (9) and (11), 

(9) H + 0 2 + M -f H0 2 + M 
(11) H0 2 + H^0H + 0H 

which contributed to the production of free radicals 
and the release of energy. This effect was very ev- 
ident for the Mach 8 condition, for which the igni- 
tion process was delayed and did not take place in 
the combustor section. Reaction still occurred in the 
nozzle section because the rate coefficient for reac- 
tion (9) increased with decreasing temperatures and 
the rate of reaction (11) was not a strong function 
of temperature. At the Mach 25 condition, omitting 
the H0 2 reactions did not affect the results. 

Thus far the discussion of the chemical kinetic 
effects has centered on the effects of groups of re- 
actions on the calculated behavior of the hydrogen- 
air system. Another aspect that needs to be con- 
sidered is the matter of the effect of the uncertainty 
in the assigned rate coefficients on the calculated re- 
sults. Baulch et al. (1972) and Baulch, Drysdale, 
and Horne (1973) have evaluated the rate coefficient 
data for most of the reactions in the hydrogen-oxygen 
mechanism and assigned error or uncertainty lim- 
its which ranged from ±20 percent to a factor of 
three. With such errors the predictive capability of 
the mechanism might be questioned. This would be 
the case if the mechanism were assembled arbitrar- 
ily without comparison or optimization with exper- 
imental data. By optimizing the mechanism with 
experimental data, we can assemble a mechanism 
with an internally consistent set of rate coefficients 
that can be defined to within rather narrow limits. 
The reaction mechanism proposed in this report was 
optimized through comparison of calculated ignition 
delay times and burning velocities with experimen- 
tal data. The resulting mechanism could reproduce 
the experimental results to within ±20 percent. As 
part of the optimization procedure, a reaction rate 
sensitivity analysis was performed in which the rate 
coefficients for all the reactions were varied within 
the error limits suggested by Baulch, Drysdale, and 
Horne (1973) and by Dixon-Lewis (1979). This anal- 
ysis revealed that the calculated ignition delay times 
and burning velocities were the most sensitive to the 
rate coefficients assigned to reactions (2) to (5), (9) 
to (15), and (6) to (8), in that order. Additional 
studies showed that varying the rate coefficients for 
reactions (2) to (5) and (9) by more than 10 percent 
about the values given by the expressions in table I 
gave calculated results that were not in good agree- 
ment with the experimental data. Similar studies 
with reactions (10) to (15) showed that variations of 
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Figure 14. Effect of HO 2 and H 2 O 2 reactions at Mach 8 condition. 
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Figure 15. Effect of HO 2 and H 2 O 2 reactions at Mach 16 condition. 




±50 percent were acceptable. Varying the rate co- 
efficients for reactions (16) to (19) by a factor up 
to two did not significantly alter the calculated re- 
sults. The most important recombination reaction, 
after reaction (9), is reaction (6). The rate coefficient 
for reaction (6) could be varied to within a factor of 
two without significantly affecting the calculated re- 
sults. The rate coefficient for the other recombina- 
tion reactions (7), (8), and (20) could be varied to 
within a factor of five and still provide satisfactory 
agreement between the calculated and experimental 
results. The sensitivity of the calculated results to 
the nitrogen-oxygen chemistry represented by reac- 
tions (21) to (33) was also examined. The calculated 
ignition delay times for the H 2 -NO-O 2 -N 2 mixtures 
(fig. 3) were very sensitive to the rate coefficient as- 
signed to reaction (30). Varying the rate coefficient 
for reaction (30) by ±50 percent did not significantly 
alter the results. The rate coefficients for the other 
reactions could be varied to at least a factor of two 
without significantly affecting the calculated results. 

Based on the results of the sensitivity study it was 
concluded that for the reaction mechanism given in 
table I, the uncertainties associated with the various 
rate coefficients are the following: for reaction (1), 
a factor of five; for reactions (2) to (5) and (9), 
±10 percent; for reaction (6), a factor of two; for 
reactions (7), (8), and (20), a factor of five; for 
reactions (10) to (15) and (30), ±50 percent; and for 
the remaining reactions, a factor of two. With these 
uncertainties only the variation in the rate coefficient 
for reaction (6), the recombination of H and OH to 
form H 2 O, would have an effect on calculated results 
for the Mach 8 and Mach 16 conditions. Varying the 
rate of reaction (6) by a factor of two produced a 
4-percent change in the chemical kinetic efficiency, 
as shown in figures 14 and 15. At the Mach 25 
condition, varying the rate of reaction (6) by a factor 
of two did not influence the calculated results. 

Concluding Remarks 

A chemical kinetic mechanism for the combus- 
tion of hydrogen has been assembled and optimized 
through comparison of the observed behavior deter- 
mined in shock-tube and flame studies with behav- 
ior predicted by the mechanism. The reactions con- 
tained in the mechanism reflect the current state of 
knowledge of the chemistry of the hydrogen-air sys- 
tem, and the assigned rate coefficients are consistent 
with accepted values reported in the literature. It 
was determined that the mechanism is capable of sat- 
isfactorily reproducing the experimental results for a 
range of conditions relevant to scramjet combustion. 

Calculations made with the reaction mechanism 
for representative scramjet combustor conditions at 


Mach 8, 16, and 25 gave the following results. Chem- 
ical kinetic effects are important and result in a de- 
crease in internal thrust compared with the equilib- 
rium chemistry condition. For Mach numbers less 
than 16, the studies suggest that an ignition source 
will most likely be required to overcome the slow ig- 
nition chemistry at these conditions. At the Mach 25 
condition, the initial temperature and pressure were 
high enough so that ignition was rapid and the pres- 
ence of an ignition source did not significantly af- 
fect reaction rates. Most of the combustion process 
occurred in the combustor section. At the Mach 8 
condition the overall reaction rate is controlled by 
the rates of the chain propagating and branching re- 
actions, the recombination reactions, and the HO 2 - 
consuming reactions. At the Mach 16 condition the 
recombination reactions and the H02-consuming re- 
actions are the rate controlling processes. At the 
Mach 25 condition the rate of the recombination re- 
action controls the combustion process. 


NASA Langley Research Center 
Hampton, VA 23665-5225 
December 22, 1987 
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Table I Hydrogen-Air Combustion Mechanism 0 


Reaction** 

A 

n 

E 

(1) H 2 + O 2 — > OH + OH 

1.70 x 10 13 

0 

48000 

(2) H + 0 2 - OH + O 

2.60 x 10 14 

0 

16800 

(3) O + H 2 — OH + H 

1.80 x 10 10 

1.00 

8900 

(4) OH + H 2 — ► H 2 0 + H 

2.20 x 10 13 

0 

5150 

(5) OH + OH — H 2 0 + O 

6.30 x 10 12 

0 

1090 

(6) H + OH + M H 2 0 + M 

2.20 X 10 22 

-2.00 

0 

(7) H + H + M — H 2 + M 

6.40 x 10 17 

-1.00 

0 

(8) H + O + M -* OH + M 

6.00 x 10 16 

- .6 

0 

(9) H + 0 2 + M — ► H0 2 + M 

2.10 X 10 15 

0 

-1000 

(10) H0 2 + H -» H 2 + 0 2 

1.30 X 10 13 

0 

0 

(11) H0 2 + H - OH + OH 

1.40 X 10 14 

0 

1080 

(12) H0 2 + H — H 2 0 + O 

1.00 X 10 13 

0 

1080 

(13) H0 2 + O - 0 2 + OH 

1.50 X 10 13 

0 

950 

(14) H0 2 + OH — ► H 2 0 + 0 2 

8.00 X 10 12 

0 

0 

(15) H0 2 + H0 2 — ► H 2 0 2 + 0 2 

2.00 X 10 12 

0 

0 

(16) H + H 2 0 2 - H 2 + H0 2 

1.40 X 10 12 

0 

3600 

(17) O + H 2 0 2 - OH + H0 2 

1.40 x 10 13 

0 

6400 

(18) OH + H 2 0 2 — H 2 0 + H0 2 

6.10 X 10 12 

0 

1430 

(19) M + H 2 0 2 -» OH + OH + M 

1.20 X 10 17 

0 

45 500 

(20) O + O + M — 0 2 +M 

6.00 X 10 17 

0 

-1800 

(21) N + N + M — N 2 +M 

2.80 X 10 17 

- .75 

0 

(22) N + 0 2 — * NO + O 

6.40 X 10 9 

1.00 

6 300 

(23) N + NO - N 2 + O 

1.60 X 10 13 

0 

0 

(24) N + OH — ► NO + H 

6.30 X 10 11 

.50 

0 

(25) H + NO + M -*• HNO + M 

5.40 x 10 15 

0 

-600 

(26) H + HNO - NO + H 2 

4.80 x 10 12 

0 

0 

(27) O + HNO — NO + OH 

5.00 xlO 11 

.50 

0 

(28) OH + HNO -* NO + H 2 0 

3.60 x 10 13 

0 

0 

(29) H0 2 + HNO — * NO + H 2 0 2 

2.00 x 10 12 

0 

0 

(30) H0 2 + NO -* N0 2 + OH 

3.40 x 10 12 

0 

-260 

(31) H + N0 2 — NO + OH 

3.50 x 10 14 

0 

1500 

(32) O + N0 2 -* NO + 0 2 

1.00 x 10 13 

0 

600 

(33) M + N0 2 - NO + O + M 

1.16 x 10 16 

0 

66000 


a The rate coefficients are given in the form k = AT n exp(-E /RT)\ units are in seconds, moles, cubic centimeters, calories, 
and kelvins. 

fc The third-body efficiencies relative to N 2 = 1.0 are as follows: for reaction (6), H 2 O = 6.0; for reaction (7), H 2 = 2.0 and 
H 2 0 = 6.0; for reaction (8), H 2 0 = 5.0; for reaction (9), H 2 = 2.0 and H 2 0 = 16.0; and for reaction (19), H 2 0 = 15.0. 


Table II. Initial Combustor Conditions 


M 

<t> 



T, K 

p, atm 

V, m/sec 

8 


670 


2090 

16 

mm 1 ■ 1 m 

■ 


4658 

25 


1 


7571 
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